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Abstract 

We calculate the self-energy of the A (1232) resonance in a finite volume, using chiral effective 
field theory with explicit spin-3/2 fields. The calculations are performed up-to-and- including fourth 
order in the small scale expansion and yield an explicit parameterization of the energy spectrum of 
the interacting pion-nucleon pair in a finite box in terms of both the quark mass and the box size L. 
It is shown that finite- volume corrections can be sizeable at small quark masses. 
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1 Introduction 



The recent surge of interest in lattice calculations of the excited baryon spectrum [1-19] has 
been mainly motivated by the experimental resonance physics program at Jefferson Lab [20] 
and ELSA [21]. Also, the hadron spectrum is arguably the least understood feature of Quantum 
Chromodynamics. In general, the extraction of the properties of the excited states from the 
lattice data is a more delicate enterprise as compared to the ground-state hadrons. The reason 
is that the excited states are unstable and, strictly speaking, can not be put in correspondence 
to a single isolated level in the discrete spectrum measured in lattice simulations. A standard 
procedure proposed by Liischer [22-25] (see also [26-30]) consists in placing the system into a 
finite cubic box of a size L and studying the response of the spectrum on the change of L. It can 
be shown that the dependence of the energy levels on L is dictated solely by the scattering phase 
shift in the infinite volume. Consequently, the method is capable of extracting the phase shift 
from the lattice data that also determines the position and the width of the resonances (see, 
e.g. [19,31,32]). Recently, the above approach has been also applied to study nucleon-nucleon 
phase shifts at low energy, as well as the two-body shallow bound states [33-37]. 

Alternative approaches to study the decaying states have been suggested, see, e.g. [38-40]. In 
particular, an interesting proposal is to reconstruct the spectral function by using the maximal 
entropy method [41], which can also be used to address the problem of unstable systems. 

In actual calculations on the lattice the quark masses do not usually coincide with physical 
quark masses. This qualitatively changes the picture since, if the quark mass is large enough, 
the A(1232) does not decay and can be extracted by the methods applicable in case of the 
stable particles. Reducing the quark mass, a value is reached such that the A starts to decay 
into a pion and a nucleon*^. The spectrum becomes strongly volume-dependent and Liischer's 
method has to be applied to extract the parameters of the resonance - the mass and the width. 

Above the threshold + M,r > the finite- volume corrections to the spectrum are 
exponentially suppressed and can be neglected in the first approximation. However, for those 
values of the quark masses which correspond to M^r -|-M^ < M^,, finite- volume corrections may 
become large and should be taken into account. Note that merely making the volume larger 
does not suffice in the case of an unstable state. Due to the potentially large corrections, the 
finite volume data on the finite- volume energy spectrum can be enhanced below threshold. This 
enhancement, which is visible in the lattice data at smaller volumes, can not be described by 
using the formulae for the quark mass dependence in the infinite volume. We shall demonstrate 
an explicit example of such a behavior below. 

From the above discussion it is clear that, in order to be able to include all available lattice 
data for large as well as small quark masses in the analysis, one needs to provide a simultaneous 
explicit parameterization of the lattice QCD spectrum in terms of both the quark mass rh 
and the box size L. This goal can be achieved by invoking the chiral effective field theory 
with explicit spin-3/2 degrees of freedom [42,43] in a finite volume. The first attempt in this 
direction was made in Ref . [44] , where we have performed the calculations of the finite- volume 
energy spectrum at third order in the so-called small scale expansion (SSE). The present paper 
extends these calculations to the fourth order. In addition, 

^"'^The decay threshold is located at Mm + Af^ = Ma in the infinite volume. In a finite volume, the decay of 
A at threshold in the center-of-mass (CM) frame is forbidden. Still, for brevity, the point + -/Vf^r = Ma will 
be always referred below to as the threshold. 
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i) We provide an explicit formula for the finite-volume corrections for the unstable A, which 
can be used in the analysis of the lattice data; 

ii) We perform a fit of the obtained expressions to the most recent available data at dif- 
ferent quark masses, taking into account finite-volume corrections. The fit allows one to 
determine some of the low-energy constants (LECs) in the chiral Lagrangian; 

iii) In doing so, one does not need to resort to any input phenomenological parameterization 
of the resonant amplitude, because SSE provides such a parameterization automatically, 
order by order in the e-expansion (here, e denotes the formal small expansion parameter 
in the SSE). 

iv) We analyze the quark mass dependence of the spectrum by using the method of probability 
distribution, introduced in [45]. 

Note also that in this paper we do not consider the finite- volume effects in the stable particle 
masses, which are exponentially suppressed at large volumes. Such effects can be treated within 
the same approach, see, e.g. Ref. [46]. 

The layout of the paper is as follows. In section 2 we discuss the calculation of the mass of 
the nucleon and the A in the infinite volume, at fourth order in the small scale expansion. In 
section 3 the calculation of the finite- volume energy spectrum of the irN system is addressed. In 
section 4 we consider the fit of the explicit analytic expressions for the nucleon and A mass to the 
existing data from lattice QCD and determine some of the LECs of the chiral Lagrangian. We 
also analyze the finite volume spectrum with the use of probability distributions [45]. Finally, 
section 5 contains our conclusions. 



The mass of the nucleon and the A resonance in the 
infinite volume 



Our calculations will be carried out in two steps. We first perform calculations of the nucleon 
and A mass at order in the infinite volume*^ . At the second step, we use the same Lagrangian 
in order to carry out the calculations of the finite- volume energy spectrum. The results of these 
calculations are applied to the case of an unstable A. 

The Lagrangian of pions, nucleons and deltas up-to-and-including order e"^ in the SSE is 
taken from Ref. [47]. Below we display only those terms that contribute to the nucleon and A 
mass at this order, 
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where the pion-nucleon Lagrangians are given by 
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ttN 



(2) 
■kN 
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rriN + — ^75 
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h.c.) + I {u') + 



N 



(3) 
■kN 



[S23Ao(x+) + + ...]^N 



"^^The small parameter e subsumes external momenta, the pion mass and the nucleon-delta mass splitting. 
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The ttA^A interaction is described by the following Lagrangians 



= CA^l.O^''w],i,N + h.C. 



C 



(2) 



TUN 



N 



h.C. . 



(4) 



In the above expressions, -j/^tv and ip^^ denote the nucleon and the A field, respectively, rhiq and 
rfiA stand for their masses in the chiral limit and Aq = rhi^ — rhiq. Note that M^^ = 0(e) and 
Aq = 0(e). The building blocks that are used in the construction of the above Lagrangian are 
given by 



D^ = d^ + -[u\d^u\ 



X = 2B{s + ip) , x± = u^X'^'^ i ux^u , s = ml + . . . 



2 
d 



(5) 



and the isospin projectors are defined by 



^3/2 



^1/2 _ J- 



(6) 



In these formulae standard notation is utilized. Namely, we use U = exp(ir ■ vr/F), where vr 
is the pion field. We work in the isospin limit m„ = = m and the trace in flavor space 
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Figure 1: Graphs contributing to the nucleon self-energy at C(e^) in SSE. Sohd, double solid 
and dashed lines denote nucleons, deltas and pions, in order. 

is denoted by (...). The quantity F is the pion decay constant, B is related to the quark 
condensate and qa is the nucleon axial- vector constant (all in the chiral limit). The coefficients 
Ci , (2^ 5 ' ' ' are the pertinent LECs. 

The propagator of a Rarita-Schwinger field in d dimensions is given by 



1 d-2 Pf/Ju^^I^ 



■ (7) 



The calculations are carried out in infrared regularization. Differently from Ref. [47,48], we 
do not project out the redundant spin-1/2 components of the A propagator, which appears in 
the loops. This amounts merely to a redefinition of some of the LECs - hence, the numerical 
values of LECs determined from fitting to the same data, should in general differ in these two 
schemes. For related discussion of this issue, see also [49-53]. 

The self-energy of the A is complex on the mass shell for those values of the pion masses, 
when the A turns unstable, i.e., M^r < Ma — Mat. The mass of the A is defined as a real part 
of the pole position in the propagator. 

The diagrams that contribute to the nucleon and A masses at order e^, are displayed in 
Fig. 1 and Fig. 2, respectively. The calculations are pretty standard and the final results are 
listed in Appendix A. Since we are primarily interested in fitting the quark mass dependence to 
lattice data, it is useful to normalize both quantities at the physical value of the quark (pion) 
mass 

= MN + x,{Ml-Ml)+X2{Ml-Ml)+x,{Mt-Mt)+xJMt\Ti^-Mt\Ti^^ 



rriN ^ M, 



N 



($A.(m;v, ^A, Ml) - $^(M^, Ma, MD) + 0{€') 



= + y,{Ml - Ml) + y^iMl - Ml) + y,{Mt - M^) + y, ( M^ In ^ - M^ In ^) 



($A(m^, mA, Ml) - <l'A(Miv, Ma, M')) + O(e^) 
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Figure 2: Graphs contributing to the self-energy of A at 0{e^) in SSE. For notation, see Fig. 1. 



where in the fit we use 

rriN = Mn + Xi{MI-MI) + --- , = + y,{M^ - M^) + ■ ■ ■ , 

2 , m'i — mlr — 

Z = + 2(mA - niMjCAh H caOq 

rriN 

= cl + 2AoCA{h + be) + -- - . (9) 

Here, = 2mB and M^r, Mat, M/^ stand for the physical values of the pion, nucleon and A 
masses. The ellipses denote the higher-order terms in e. Further, at the order we are working, 
one may take Aq = Ma — + ■ ■ ■ in the above equations. The masses Mat, Ma are functions 
of the pion mass M^. The quantities Xi,yi, Z denote certain combinations of LECs. Explicit 
expressions for the Xj, yi, Z, as well as for the functions $Ar, $a are displayed in Appendix A. 
Fitting the nucleon and A masses, given by Eq. (8), to the lattice data determines the numerical 
values of the above combinations of LECs. Note that some higher-order terms are also present 
in Eq. (9), e.g. in the expressions for $a, '^'Af- 

The calculation of the quark mass dependence of the nucleon and A masses has been carried out 
in different settings [47,48,54-56]. Note that, in particular, our result for the nucleon mass in 
the infinite volume agrees at 0{e^) with the expression given in Eq. (17) of Ref. [54]. However, 
it differs from the O(e^) result for the nucleon and A masses, which are displayed in Eqs. (22) 
and (30) of Ref. [55], respectively. For instance, these latter expressions do not contain the 
LECs which describe the quark mass dependence of the uNA vertex (analog of the constants 



3 Self-energy of the A resonance in a finite volume: the 
energy levels 

3.1 Calculation of the finite-volume correction 

In a finite volume the A propagator develops a tower of poles on the real axis. The location of 
these poles determines the finite- volume energy spectrum of the system. Thus, calculating the 
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propagator in a finite volume, we shall be able to study the volume- dependence of the energy 
levels. The procedure is described in detail in Ref. [44] and will not be repeated here. Here 
we simply note that the only difference to the infinite- volume case is the replacement of the 
(Euclidean) loop integrations by infinite sums 

In the above expression, L denotes the size of the (cubic) box in which the system is placed. 
The Lagrangian that produces the loops is the same as in the infinite volume. 

The calculations are substantially simplified, if carried out in a large volume where the 
exponentially suppressed corrections can be neglected. In this limit the masses of the stable 
particles can be considered as volume-independent. However, as it is well known, the energy 
levels corresponding to the unstable particles receive corrections, which are suppressed by pow- 
ers of L. Only the diagrams, which contain the pion-nucleon intermediate state - diagrams 
b,e,f,g in Fig. 2, calculated in a finite volume - contribute to this power-like behavior. Retain- 
ing the finite-volume parts of these diagrams only, the equation that determines the location 
of the poles in the A propagator is written as (cf with Ref [44]) 

M^-E=^^{iE + M^f - Ml) W,^{E% . (11) 

Here, E denotes the pole position on the real axis, and M^, M^v, Ma are the masses in the 
infinite volume. Further, Z stands for the following combination of the LECs 

Z = c\ + 2hcA{E-MN) + 2hcA — f ^ 

= Z + 2{b, + b,)cA{E-M^) + 0{e'). (12) 

It is seen that only three LECs: Ca,&3,&6 appear in the finite-volume correction to the energy 
of A. 

Finally, the quantity W^lE"^) corresponds to the finite- volume part of the nN loop function 



W^' {E') = W^' {E') - ) , (13) 



where 



W^^(E') = f^l-y 1 ^ , A = {^E, 0) . (14) 

°^ W 27rL3 Y(M2 + F)(M2, + (P-A;)2)' a« ^ ' ^ ^ ' 

In large volumes, neglecting exponentially suppressed contributions, the loop function above 
threshold can be rewritten as 



Zoo(l,g^) = 2:00(1, g^) - 2^00(1, g^' 



(15) 
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Figure 3: Feynman diagram yielding the same scattering phase as Eq. (18) from an infinite 
volume SSE calculation. Solid, dashed and double lines denote nucleons, pions and deltas, 
respectively. 



where the ellipses stand for the exponentially suppressed contributions, the quantity q = 
with p = \^/'^{E'^, Mlf, Ml) /2E and Zqq is the zeta-function from Ref. [23] 

Note that Zoo(l, q^) = 2oo(l, q^) for q^ > 0. 



3.2 Relation to Liischer's formula 

By using Eqs. (15) and (16) it can be checked that the Eq. (11) which determines the position 
of the pole in the propagator, can be rewritten in the form of Liischer's equation 

tan5(p) = — — — , (17) 
2oo(l,r) 

where 6{p) denotes the scattering phase shift in the Pas-channel for the following choice of the 
scattering phase 

which corresponds to the s-channel tree-level scattering amplitude in the SSE, shown in Fig. 3. 
The discrete solutions of Eq. (17) determine the energy spectrum of the system = \/ Mjj + p^ 
+ a/ -|- pI^ through the given scattering phase S{p). 



3.3 Effect due to the finite lattice spacing 

Certain caution is needed, if one uses the above formulae in order to fit the lattice data. Indeed, 
they contain artefacts due to the finite lattice spacing a. For example, in the analysis of the 
data obtained by using twisted mass fermions, one has to address the issue of isospin breaking 
at finite a. Even if the effect turns out to be not very large in the measured nucleon and delta 
masses, the neutral pion masses in the loops will differ strongly from the charged ones. It is 
clear that, in order to address the problem in its full generality, one has to develop twisted 
mass chiral perturbation theory, where the isospin breaking emerges at a finite lattice spacing. 
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In this paper, however, we shall restrict ourselves to the spectrum of A'^"'",A~, where only 
charged pions occur in the loops up-to-and-including order e^. Consequently, at this order one 
may use the conventional formalism, with the pion mass set equal to the charged pion mass and 
assuming isospin symmetry in the couplings. The data on A"*", A° will be used for checking the 
size of isospin-breaking contributions at finite a and thus will serve as an error estimate only. 



3.4 Determination of the width 

The width of the A at the physical value of the quark mass is determined by the parameter Z 
which, in turn, at this order depends on the LECs c^, 63, h^, see Eq. (9) 

Zql^ {M^ + Mm?-MI 

= mi ■ ^^^^ 

In the above equation, Qcm denotes the CM momentum of the vrA^ pair after the decay of A, 
is the pion decay constant and Ma = 1232 MeV. 

A determination of the LECs c^, &3, from the fit to the A mass in the infinite volume 
does not provide sufficient accuracy, because these LECs enter starting from the next-to-leading 
order. The situation changes, however, if we consider the data obtained at the same quark 
mass and at different volumes. Consider, for instance, the data taken at two different values 
of L. Since the mass of the A in the infinite volume is, by definition, volume-independent, the 
following consistency condition must hold at this order 

Ma = ^a(^i) + 5E^{L^, CA, h, h) = E^iL^) + 5^a(^2, ca, ^3, h) , (20) 

where E^iLi), i = 1,2, denote the measured energies and 6E/^{Li, ca, bs, b^) denotes the finite- 
volume correction below threshold evaluated at the pertinent values of E and L, see Eq. (11). 
Performing measurements at different values of L provides additional constraints. Extracting 
the values of LECs from the above conditions, one may in principle determine the width of the 
A by using Eq. (19). Note that Eq. (20) holds at a fixed value of the quark mass. 



4 Fit to the lattice data 
4.1 Choice of the data 

Just in order to demonstrate the application of the theoretical framework developed above, we 
shall perform the fit to the recent data of the ETM collaboration [12]. In particular, we fit the 
data for the nucleon and A masses, obtained on (3 = 3.9 lattices of size 24^ x 48 and 32'^ x 64 
(smeared link and smeared source), corresponding to L = 2.1 fm and L = 2.7 fm, respectively. 
These data are given in Table II of Ref. [12]. The data contain the nucleon and A masses at 4 
different values of the quark mass (on a smaller lattice) and one additional data point for the 
lightest quark mass (on a larger lattice). At the lightest quark mass, the sum of the nucleon 
and pion masses is smaller that the A mass. Note that we do not have access to the data at 
different volumes, extrapolated to the continuum limit a — >■ 0. The values of the nucleon and 
A masses, displayed in Table II of Ref. [12] still contain the artefacts due to a finite lattice 
spacing. 
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4.2 Fit to the nucleon and A masses: infinite volume 

In Ref. [12] the infinite-volume mass of tlie A is identified with the extracted energy level at 
a largest volume at a given quark mass. As already mentioned, such a procedure can not 
be strictly justified for unstable particles. Notwithstanding, we shall use this method in the 
beginning and try to simultaneously fit both nucleon and A masses with the infinite-volume 
formulae (8). The result is shown in Fig. 4. For comparison, in the same figure we display the 
data points taken on a smaller lattice. 

The Eqs. (8) contain too many free LECs, making the fit to the few available data points 
questionable. A reasonable strategy consists in constraining some of these LECs by using 
additional physical information. Thus, the LEG X2 is unambiguously fixed through the known 
value of the nucleon axial-vector coupling qa = 1.267. Furthermore, we use the S't/(6)-relation 
Qi = (9/5)(7a and set 02,3 = 02,3. The LECs 02,3 are determined by matching to ChPT without 
explicit A degree of freedom. The pertinent relations are given by C2 = C2— g^/ {2Ao) +0{1) and 
C3 = £3 + ^y(2Ao)+0(l), where 02,3 denote the LECs in ChPT without A (cf with Ref. [53]). 
Using the values C2 = 3.3 GeV"^ and C3 = -4.7 GeV^^ [52], we finally get C2 ~ 0.55 GeV"^ 
and C3 ^ —1.95 GeV~^. In addition, we use the value Z = 2.14 that leads to the physical decay 
width r = 118 MeV after substituting into Eq. (19). The couplings ca and 63 -|- Bq are given 
below, see Eq. (22). 

The remaining LECs Ci,ei,ai,ef are allowed to vary freely (these LECs are defined in 
Eq. (A.3)). In the fit we will use the data at L = 2.1 fm except the lowest point corresponding 
to L = 2.7 fm. The fit to the data gives the following values for these parameters (no errors 
assigned) 

ci = -1.6 GeV"^ , ai = -1.8 GeV"^ , 

ei = -1.1 GeV~^ , ef = 6.6 GeV"^ . (21) 

The SU{6) relation di ~ ci holds approximately, in difference with the result obtained in 
Ref. [47] (note, however, that the different prescriptions for performing the infrared regulariza- 
tion in the case of A amount to a finite renormalization of various LECs). In order to compare 
the obtained value of ci with the phenomenological estimates, one has again to perform the 
matching to ChPT without an explicit A, which yields Ci = Ci + ZAo/(87r^F^) ln(2Ao/m7v) + 
0(Aq). In this expression, ci denotes the value of the pertinent LEG in ChPT. The resulting 
shift ~ 0.44 GeV~^ in ci is positive, and the obtained value of Ci reasonably agrees with the 
value extracted from the phenomenological analysis of the pion-nucleon scattering at fourth 
order, see e.g. [57] for the latest update. Note, however, that the fourth order LECs ei and 
ef differ significantly. Moreover, ef is rather large that could serve an indication of a poor 
convergence at these pion masses. 

As one observes from Fig 4, the finite-size corrections to the A energy may turn out sizable 
below threshold (the data point corresponding to the smallest pion mass). At present, the error 
bars on the data are large that precludes one to make an unambiguous statement on the issue. 
However, even at the present accuracy a hint is seen that the lowest data point at L = 2.1 fm 
is located above the curve. This is an example of the enhancement which was mentioned in the 
introduction. 
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Figure 4: The fit to the nucleon and A"*"^ spectrum by using Eq. (8). The filled circles cor- 
respond to the data taken at L = 2.1 fm. The data corresponding to L = 2.7 fm at the 
smallest pion mass are shown for comparison (triangles). The black diamonds without error 
bars correspond to the physical masses. 

4.3 Analytic behavior at threshold 

It is quite instructive to study the qualitative behavior of the energy levels in the vicinity of 
threshold, i.e. choosing the quark mass so that the sum of the pion and nucleon masses are 
only slightly below the A mass. As we know, this situation is realized for the lowest data point. 

Let us consider the plot of the function g^^oo(l; (f)-, which enters the r.h.s. of Eq. (11), see 
Fig. 5. This quantity has a cusp, proportional to g^, at threshold = 0. Moreover, its value in 
the limit — > is different from zero. Below threshold, the function decreases exponentially. 
Above threshold, the function has a tower of poles, with the first one located at = 1. 

If one is varying the quark mass so that g^ stays negative (A stable), the finite- volume 
corrections are exponentially small. However, if decreasing the quark mass, the quantity g^ 
moves across the cusp from below, the effect blows up rapidly. In this case, the energy levels 
in a finite volume receive large corrections, which should be taken into account. On the other 
hand, the "raw" data on the energy levels at a fixed volume, which are depicted, e.g. in Fig. 4, 
are smooth functions of the quark mass and do not exhibit any cusp. 

4.4 Subtracting finite-volume effect 

In order to subtract finite-volume effect at order e^, one has to fix the values of the LECs 
and 63 + fee- Since we have only one data point below threshold, both LECs can not be fixed 
simultaneously. For this reason, we have set the constant Z = 2.14 so as to reproduce the width 
of the A and used the consistency condition (20) to determine Z and thus to disentangle ca 
and 63 + 65 from Eqs. (12) and (A. 4). Using central values for the energy levels, we get 

4 = 2.73 , 63 + fee = -0.6 GeV"^ . (22) 
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Figure 5: The function q"^ Zqq{1] q"^) in the vicinity of threshold = 0. The threshold cusp is 
clearly visible. 



L [fm] 




Mm 


E^++.- 








2.1 


314 ±2.4 


1189 ± 14 


1574 ± 29 


1609 ± 40 


-90 


-129 


2.7 


309 ±1.9 


1177 ± 13 


1523 ± 23 


1523 ± 34 


-39 


-43 



Table 1: Meson and baryon masses for two different values of the box size L (the data are 
taken from Table II of Ref. [12] (central values only) and correspond to the choice SS of the 
interpolating field). Last two columns correspond to the finite- volume corrections to the energy 
levels, calculated by using Eq. (11) (see the text for more detail). All masses are given in MeV. 

As seen, these LECs are indeed of the natural size. 

In Table 1 we give the results for the finite-volume corrections to the central value of the 
lowest data point, evaluated at the above values of the LECs. These finite volume corrections 
are indeed small except for the lowest point. The results for A"'"'° are presented just for the 
visualization of the artefacts due to the finite lattice size. As is seen from this table, the finite- 
volume corrections matter even at the present accuracy. For instance, the infinite-volume mass 
of the A"'"'' is equal to 1484 MeV. Here we note that in Ref. [58] significant finite-volume 
corrections have been found as well. The calculations in Ref. [58] have been carried out at 
order e^, by using the formula of Ref. [44]. At this order, one would set = Z = 2.14 and 
^3 ± ^6 = in our formulae. It can be checked that this does not change the result significantly. 

In Fig. 6 we show the fit to the lattice data. The finite- volume effect, which is given in 
Table 1, is subtracted from the lowest data point. It can be seen that the LECs, which are 
extracted from the fit, are quite stable (to be compared to Eq. (21)) 

ci = -1.6 GeV"^ , ai = -1.7 GeV"^ , 

ei = -1.4 GeV"^ , ef = 6.4 GeV"^ , (23) 

however, is somewhat worse in this case. 

As can be observed from Fig. 6, the finite-volume correction to the lowest data point is 
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Figure 6: The fit to the nucleoli and A++ spectrum. The lowest data point for A has been 
purified with respect to the finite-volume corrections. For comparison, the uncorrected lowest 
data points for L = 2.1 fm and L = 2.7 fm (triangles) are shown. 

significant. There is no enhancement in the corrected data. 

Finally, just as a hint, we would like to mention that it is possible to get a very good fit to 
the data, concentrating only on two lowest quark mass data points and relaxing the condition 
gi = {9/5)gA- The obtained values for the LECs are (71 = 2.89 ~ 2.3gA, ci = —1.43 GeV~^,ai = 
— 1.67 GeV~^,ei = —1.35 GeV~^ and ef = 2.02 GeV"'^. As can be seen, ef is now of natural 
size. The constant ci, contributing to the nucleon a-term at lowest order, turns out to be 
slightly smaller. Of course, two data points do not provide sufficient input to draw definite 
conclusions about the values of the LECs. For the same reason, we refrain here from citing the 
values of the nucleon and delta a-terms, which can be reliably determined, only if more data 
points become available at smaller quark masses. 

4.5 Probability distribution: dependence on the quark mass 

In this section we shall study the quark (pion) mass dependence of the structure of the energy 
levels. To this end, it is useful to invoke the language of the probability distributions [45], 
which makes this dependence very transparent. 

The probability distribution, which can be constructed from the volume-dependent energy 
spectrum through an unambiguous procedure [45], is closely related to the so-called density 
of states in a finite volume. Using Liischer's formula, it can be shown [45] that - to a good 
approximation - the probability distribution W{p) can be expressed via the scattering phase 

W(p) = - V (VE^^HM + 27r^-(p) 

P tl\ P VM^n - 6ip)) 

where 6{p) denotes the scattering phase, is the number of energy levels analyzed and C 
denotes the normalization constant. Below we restrict ourselves to the analysis of the lowest 
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p, GeV 

Figure 7: Subtracted probability distributions for different values of the pion mass. The quan- 
tity p is the relative momentum of the ttN pair in the center-of-mass frame. The solid lines 
correspond to the theoretical prediction based on Liischer's formula, see Eq. (24). 

state, putting = 1. 

In case of the wide resonance like A, it is convenient to consider the so-called subtracted 
probability distribution, which is obtained from W{p) by subtracting the background Wfree{p) 
corresponding to the free vrA^ pairs with 6{p) =0 [45]. In the vicinity of the resonance, the sub- 
tracted distribution approximately follows the Breit-Wigner form of the scattering cross section 
and thus allows one to easily identify the resonance from the data on the energy spectrum. 

Using the values of the various LECs determined from the fit, and substituting the scattering 
phase given by Eq. (18) into Eq. (24), one may easily predict the shape of the probability 
distributions at different values of the pion mass. The results are given in Fig. 7. It is seen 
that the distributions behave in the expected manner: for the higher pion masses, the center- 
of-mass momentum decreases and the distribution becomes narrower. Slightly after 300 MeV 
the distribution degenerates into the (5-function - the A resonance becomes stable. Of course, 
such a behavior can be only observed in practice, provided there are at least few data points 
with different values of L at a given quark mass [45] . 

5 Conclusions 

Here, we summarize the pertinent results of our study: 

i) We have calculated the ground state energy of the A resonance in a finite volume up-to- 
and- including O(e^) in the small scale expansion. 

ii) The obtained explicit expressions have been used to analyze the recent data on the nucleon 
and A spectrum, provided by the ETM collaboration [12]. It turns out that the finite 
volume corrections are sizable using the central value for the data point with the smallest 



14 



quark mass. Even at the present accuracy, this correction should be taken into account. 

iii) It is checked that the numerical values for the correction at O(e^) and at O(e^) do not 
differ significantly. 

iv) We perform a simultaneous fit to the nucleon and A masses in the infinite volume. The 
values of the LECs obtained in a result of such a fit are stable. However, the convergence 
of the chiral expansion of the baryon masses in the infinite volume is rather poor, since 
we are still at relatively high values of the quark masses. 

v) The measurement of the energy spectrum at different volumes opens the possibility for 
the extraction of the decay width. To this end, we have proposed a procedure based on 
the consistency condition Eq. (20). This procedure can be used meaningfully, provided 
more accurate data points emerge below threshold. 
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A The masses of the nucleon and the A(1232) 

A straightforward calculation of the nucleon and A masses at fourth order yields 



3M4(8ci-C2-4c3) M , f „ , ^„ , M 

In h 5- < Pi + 2P2 In 



32'k'^F'^ rriN ASir'^F'^mNm'^ \ tjin 

((mA + m^)2-M2)A(mi,774^if2) 



QmNF'^m\ 



(mA - m^) (mA + m^)3 / (m^ - mlf M 



M'^{2m\ + 2m^ — 171/^171^] 



3m% 



O(e^), (A.l) 
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bgW 5glM^ /20 , M 49\ 3MV 
Ma = niA- ^„ - ^ In ^ TT ' 



967rF2 1927r2F2mA \9 mjv 54 y 1287r2F2 

3M^(8ai - as - 403) , M ^, M 

^ 2 In — + 5 1 Qi + In 



24m^F2 



2M2(2m^ + 2m^ - mAmN) 



3m\ 



+ 0{e'), (A.2) 



where = 2i?m and the "tree-level masses" are given by 

niN = TTiN- 4:CiM^ - 4S23A0M2 - 332^1 - EiA^ - AE^AlM^ - Ae^M^ 

= m^v - 4ciM2 - 4eiM\ 
mA = niA- 4aiM^ - 45f AqM^ - B^^l - E^A^ - AE^AIM^ - Ae^M^ 

= mA - 4aiM2 - 4ef M\ (A.3) 
where ei = 4e38 + | (ens + ene) and ef = 4e^ + ^ (ef^g + ef^g). Furthermore, 

2 m'i — mi, — 
Z = + 2(mA - mN)cAh H c^Oe 

= 4 + 2Aoc^(63 + &6) + 0(e2), (A.4) 



and 



Pi = — <^ ^ -3m2, - 3mAm^ + 8m^ 

2mlr 3 



4,0 ^ 2 , 2 \\ I I, , '"^N + "^A 



- (3m^ + 2mA(mA + m^)(mA-mAr))cA &3+ „ 
P2 = c\{-m\ + m\ + 3mAmAr) 

\ / 2 2 f -1 f^N + TTlA 

- (mA + mN)(3m^ + 2mA-2mArnN)cA{h-\ 

\ IrriN- 

Qi = c^(3mA + Am% + AmlfniA + Am\m\ + AmMm\) 

I , \fo h , i^A + mN)cAbG 

- (mjv + mAj 2c^63 + 



rriN 

X (3mA + 4mAr(m^ + mA)(mAr - rriA)) , 
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at; 



(rriN + niA) 2cAh + 



(mA + mN)cAbe 

TUN 



-2mNmA + Sm\ + 2m^) 



The loop functions are given by 



-X 



167r^r 



^ arccos 

'JV 
,2 I A,f2 



327r2m^ 



21n^ 



if A < 



In 



21n^ 



if A > 



W, 



N 



m, 



■ -.y^^i arccos f - 



327r2m^ 



21n^ 



if A < 



' 327r2m^ 



2^ +M2-m2,-VA 



'A '"-iV 
327r2m^ 



21n^ 

mjv 



if A > 



where A = X{m\, m%, M^). 

We further express the quantity through the pion mass, according to 



Mi 1 + 



Mi 



327r2F2 



/3 + In 



Ml 

M2 



(A.5) 



(A.6) 



(A.7) 



where I3 = 2.9 ± 2.4 is the O(p^) LEG in the meson sector of chiral perturbation theory and 
stands for the physical pion mass. 

Finally, normalizing Mn and Ma at = and neglecting higher-order terms in the 
e-expansion, we obtain the equations (8) from section 2, where 



Xi 

yi 
2/2 



-4ci, 
-4ai , 

327rF2 ' 

967rF2 ' 
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= _4 A _ 49 3a2 

1927r2F2mA 54 1287r2F2 



t2 



I'M 



3(74 3(8ci — C2 — 4C3) 1 , ., 

^4 = - QO 2F2- + 00 2F2 + g 2 F2 - " 2CAmjV &3 + 06 



2/4 



47r2F2 

5gl 20 3(8ai - aa - 4a3 



1927r2F2mA 9 327r2F2 

1 fr^ 2 ^c^ - /I . 1 \\ . '^1 



7c^ - UcArriAibs + feg 



1927r2F2mA V ^ ^ ' oyy . ^^2^2 ' 

and the "tree-level" masses mM,m/\ are given by Eq. (9). 
Finally, the loop functions in Eq. (8) are defined as 

^ I ..2, ((mA + - Ml)\{ml, m%, Ml) ^^^^ ^ 



+ 



N 



Ml {2ml + 2^2, - mAm^v) 



3m|r rriM 
^ ^ - (("^A + niNf - Ml)X{m\, ml, Ml) , n 



+ 



(mA - mM){mM + m/^)'^ ( {m\ - m^)2 



3847r2m^ \ 3m^ 

2M2(2mi + 2m^-mAm^) ^ ^^p^^M^ 
3m^ mN 
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